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Abstract. Percolation is perhaps the simplest example of a process exhibiting a phase 
transition and one of the most studied phenomena in statistical physics. The percolation 
transition is continuous if sites/bonds are occupied independently with the same probability. 
However, alternative rules for the occupation of sites/bonds might affect the order of the 
transition. A recent set of rules proposed by Achlioptas et al. [Science 323, 1453 (2009)], 
characterized by competitive link addition, was claimed to lead to a discontinuous connectedness 
transition, named "explosive percolation". In this work we survey a numerical study of the 
explosive percolation transition on various types of graphs, from lattices to scale-free networks, 
and show the consistency of these results with recent analytical work showing that the transition 
is actually continuous. 



1. Introduction 

Phase transitions lie at the heart of the modern development of statistical physics They 
are changes in the state of order of a system and can be classified based on their properties in the 
immediate proximity of the critical point. Phase transitions are continuous (or second-order), 
if the order parameter changes continuously across the two phases, with an infinite correlation 
length and consequent power law decay of correlations and divergence of higher moments of the 
order parameter at the critical point. Otherwise one speaks of discontinuous (first-order) phase 
transitions, which are typically characterized by a discontinuous jump of the order parameter 
at the critical point. 

The purely geometric process known as random percolation [3\ offers a paradigmatic example 
of a continuous phase transition. The starting point is a graph, e. g. a lattice. The sites or 
the links of the lattice are occupied independently with some probability p. Nearest-neighboring 
occupied sites/links form structures called clusters. For low values of the occupation probability 
p just a few small clusters are formed, but if p increases the number and size of the clusters 
will increase as well. When the occupation probability exceeds a critical value pc, a macroscopic 
cluster, occupying a finite fraction of the total number of sites/links, emerges. This macroscopic 
structure is called percolation cluster and its relative size P, the percolation strength, is the order 
parameter of the transition: P = indicates the phase with only microscopic clusters, whereas 
P > indicates the phase with (at least) one macroscopic cluster. Random percolation has 



been studied on lattices, random graphs [1] and scale- free networks [3 El [71 [H [9] . Analytical and 
numerical studies have proved that the percolation transition is continuous, without exceptions. 
This however holds for random percolation. It cannot be excluded a priori that alternative 
processes of occupation of sites/links might lead to different types of geometric transitions. 

In a recent paper |10] . Achlioptas et al. have introduced a special set of rules, in which links 
are occupied as a result of a competitive process between pairs of links {Achlioptas processes). 
The idea is to slow down the process of cluster growth, by inserting links leading to the merge 
of small clusters. This can be done in several ways. Achlioptas et al. focused on the so-called 
product rule (PR): given a pair of links, randomly selected among those which are not yet 
occupied, one occupies the link merging the two clusters with smaller product size (Fig. [T]). 
One could consider variants of this rule, like taking the sum instead of the product, or just 
the minimum size of the clusters of each pair. Also, the competition can be extended to more 
than two links. In any case, the result of such processes is a slow growth of the cluster sizes, 
which causes a delay in the onset of the percolation transition. On the other hand, since the 
density of links and, consequently, of clusters at the onset is higher than for random percolation, 
it is natural to expect that the percolation cluster has a very rapid growth. This is indeed 
confirmed by numerical studies; in fact, the growth of the percolation cluster is so quick that 
the percolation strength P appears to vary discontinuously at the onset. The sudden jump in 
the order parameter has motivated the name "explosive percolation". In the last two years, 
Achlioptas processes have been extensively studied and meanwhile a lot is known about the 
explosive percolation transition [IIlliainilllliaiiailZllISlliaEQlEIlES]. 

The main issue concerned the order of the transition: is it continuous or discontinuous? 
Achlioptas et al. claimed that it is discontinuous, and this is what has mostly attracted the 
attention of scholars. As a matter of fact, it was soon shown that, despite the alleged jump 
of the order parameter at the threshold pc, the explosive percolation transition has peculiar 
features of continuous phase transitions, like power law distributions of cluster sizes [151 El HH] 
and power law scaling of the mean cluster size at pc [121 [131 [IS]- Indeed, in recent works by 
Nagler et al. pT] and da Costa et al. [22] it was proven that the transition is actually continuous, 
and that the jump of the order parameter is only apparent, and due to the very small critical 
exponent of the order parameter /3. 

In this paper we summarize the finite size scaling analysis of Ref. [15], and verify that the 
results are indeed consistent with the analytical findings of da Costa et al.. In Section [2] we 
introduce finite size scaling, in Section [3] we discuss the results for different types of graphs. A 
summary is reported in Section [H 



2. Finite size scaling 

Finite size scaling [23J is a well-known technique adopted in numerical studies of phase 
transitions. For continuous phase transitions, the correlation length is infinite at the critical 
threshold pc, so every variable X is scale-independent in the proximity of that point and has a 
power law form, 

x^\p-Pcr, (1) 

where a; is a critical exponent. If the system has a finite size A'^, the variable X near the threshold 
has the following scaling form 

X = N-'^/^F Up - Pc) A^^/'^j . (2) 

In Eq. [21 is a critical exponent and F a universal function. Exactly at the critical point 
ijp = Pc) the variable displays the simple scaling X ~ which can be used to extract the 

exponents' ratio uj/u, by using systems of different sizes. Moreover, if and u are known, 







Figure 1. Scheme of an Achlioptas process with product rule. Two links (dashed lines) are 
selected at random among all possible pairs of non-adjacent nodes. The link joining the pair of 
clusters with the smaller product size is eventually added to the system. Here the winning link 
is that between clusters ci and C2 (product size 7-2 = 14<4-4 = 16). 



the expression XN^I^ as a function of (p — pc) N^^'^ is just the universal function F ^ which is 
independent of A'^, so curves referring to different system sizes collapse. 

We investigated the two main variables of percolation [3], i.e. the 'percolation strength P 
and the average cluster size S. The percolation strength P, as we have said above, is the order 
parameter of the transition, and equals the relative size of the percolating cluster(s) with respect 
to the total system size N. While on lattices there are operative criteria to define a percolating 
cluster (e.g. if it runs from one edge of the lattice to the opposite one, say), on generic graphs 
this is not the case, so P is defined as the relative size of the largest connected cluster. The 
scaling ansatz of the percolation strength is 



The average cluster size S is defined as 



S 



(3) 



(4) 



where ris is the number of clusters of size s per node. The sums run over all possible values of 
s except for the one of the largest cluster. The scaling ansatz of S is 



S = Ar7/^i7(2) Up _ ^l/u 



(5) 



The universal functions F^^'^ and F^'^^ of Eqs. ^ and ([5]) are not the same, but they are related. 

In random percolation, the probability distribution P{s) of cluster sizes (except the largest), 
decays at pc as the power law P{s) ~ s~'^ with the cluster size s. We have computed the cluster 
size distribution P{s) at pc and measured the Fisher exponent r. For a given system, P{s) is 



related to Ug by the relation P{s) = Nns/uc, where ric is the total number of "finite" clusters. 
We shall use the symbol Ug to indicate P{s) as well, but in the plots Ug is normalized as P{s), 
for consistency. 

The percolation threshold pc is localized in two independent ways. The first method exploits 
the scaling of the pseudo-critical points Pc{N) 

Pc = Vc{N) + 6iV-i/- . (6) 

By using several system sizes, one can perform a fit with the three parameters 6, v and pc- The 
pseudocritical point for a system with finite size can be defined in several ways, for us it 
indicates the value of p at which the average cluster size S peaks. 

An alternative procedure relies on Eq. ([3]). The percolation strength P is plotted as a function 
of the system size for a given value of p. Since, for p = Pc the scaling follows a power law, 
the correct value of the percolation threshold can be determined by finding the value of p which 
yields the best power law fit. 

3. Results 

3.1. Implementing Achlioptas processes with product rule 

The starting point is just a graph with nodes and no links. Links are added one by one, 
according to the competitive rule previously described, i.e. by selecting each time a pair of 
links at random and picking the one yielding the smaller product for the sizes of the clusters 
it merges. For scale-free networks the situation is a bit more involved, and we describe it in 
Section 13.41 The procedure goes on until the desired density of links p is reached. We defined p 
as the number of links of the graph divided by the total number of links present in the graph 
when it has been "completed", i.e., after the addition of the last link. All graphs considered here 
are "sparse", i.e., the average degree (k), expressing the ratio between (twice) the number of 
links and the number of nodes A^, does not depend on A^. At time t of the growth process there 
are exactly t links in the system: their density p, according to our definition, is then t/{N{k)). 

3.2. Lattices 

Achlioptas processes on the square lattice were first studied by Ziff [11], who found similar 
properties for the explosive transition as Achlioptas et al. had found for Erdos-Renyi random 
graphs. The criterion to assess the nature of the transition was the same one proposed 
by Achlioptas et al., namely the scaling with the system size A^ of the transition window 
Ap = p2 — Pi, where p2 is the lowest value of p for which P > 0.5 and pi the lowest value 
of p for which P > 1/y/N. In a recent paper Ziff has performed a finite size scaling analysis 
as well [18| . The results of our analysis are shown in Fig. [2j The order parameter P yields a 
very small exponents' ratio P/v, compatible with zero [0.07(3)] (Fig. [21 left), which is what one 
would expect to find for a discontinuous transition. The average cluster size S*, instead, has a 
non-trivial power law scaling at pc, with exponent 7/1^ = 1.7(1) (Fig. [2l center). This had first 
been observed by Cho et al. in scale-free networks [12] . Fig. [2l(right), showing the distribution 
of sizes Us for all clusters except the largest one, provides an explanation of the scaling of S. 
The distribution is a clear power law [exponent 1.9(1)], which is incompatible with a classic 
discontinuous transition, as it usually occurs for continuous transitions. Since Ug is a power law, 
all variables derived from n^, including the average cluster size S, have power law scaling. 

In 3d-lattices, the general picture is consistent with that in two dimensions (Fig. [3D. Again, 
the scaling at pc of the order parameter P yields a very small exponent P/v = 0.02(2), compatible 
with zero (Fig. [S] left). Still, S scales with an exponent ^ /v = 2.1(1) (Fig. [S] center), again due 
to the power law shape of the distribution of cluster sizes (Fig. [3] right). The Fisher exponent 
T = 1.99(4) is compatible with the value we measured in two dimensions [1.9(1)]. 



□ □ □ □ 

■(T - O- -O- _ o 



6 7 
ln(L) 



12 
^ 8 



4 5" 



.0 
✓ ♦ 



6 7 
ln(L) 



10 



10" 



10"' 



10 



10"' 



10" 



10^ 



10^ 



Figure 2. Explosive percolation on 2d-lattices. (Left) Percolation strength P as a function 
of the lattice side L for three different values of the occupation probability: p = 0.5256 (violet 
diamonds), p = 0.5266 (orange circles) and p = 0.5276 (grey squares). The dashed line indicates 
the best fit obtained at the critical point p = Pc = 0.5266(2), from which we get P/v = 0.07(3). 
(Center) The average cluster size S as a function of the lattice side L for the same values of p used 
for P. The dashed line has slope ^/f = 1-7(1). (Right) Cluster size distributions measured at pc 
for PR (grey squares) and random percolation (orange circles). The exponents are r = 1.9(1) for 
PR (red dotted line), while for random percolation r = 2.05(1) (black dashed line). Simulations 
have been performed on lattices with side L = 4096. 
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Figure 3. Explosive percolation on 3d-lattices. (Left) Percolation strength P as a function 
of the lattice side L for three different values of the occupation probability: p = 0.3866 (violet 
diamonds), p = 0.3876 (orange circles) and p = 0.3886 (grey squares). The dashed line indicates 
the best fit obtained at the critical point p = Pc = 0.3876(2), from which we get fi/v = 0.02(2). 
(Center) The average cluster size 5* as a function of the lattice side L for the same values of p 
used for P. The dashed line has slope ^/v = 2.1(1). (Right) Cluster size distributions at the 
critical point for PR (grey squares) and random percolation (orange circles). The exponents 
of the power law fits are r = 1.99(4) for PR (red dotted line) and r = 2.20(1) for random 
percolation (black dashed line). Simulations have been performed on lattices with side L = 256. 



3.3. Erdds-Renyi graphs 

Fig. m presents the results of our analysis. For PR we see again a flat profile of the order 
parameter P with (/3/i^ = 0.02(1)), consistent with a discontinuous transition, along with a 
power law scaling of the average cluster size S, with exponent 7/1/ = 0.48(4). The exponent 
T = 2.08(5) (Fig. m right) is still compatible with the values found for PR on lattices, both in 
two and three dimensions (see Section [3.21 and Table 1). 

3.4. Scale-free networks 

Many natural, social and man-made systems, if represented as graphs, display common features. 
The most striking is a broad distribution of the degree of the nodes, with a large majority of 
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Figure 4. Explosive percolation on Erdos-Renyi random graphs. (Left) Percolation strength 
P as a function of the network size N for three different values of the occupation probability: 
p = 0.8872 (violet diamonds), p = 0.8882 (orange circles) and p = 0.8892 (grey squares). The 
dashed line indicates the best fit obtained at the critical point p = Pc = 0.8882(2), from which 
we get P/v = 0.02(1). (Center) The average cluster size S as a function of the network size N for 
the same values of p used for P. The dashed line has slope ^jv = 0.48(4). (Right) Cluster size 
distributions measured at the critical point for PR (grey squares) and random percolation 
(orange circles). The exponents are r = 2.08(5) (red dotted line) for PR and r = 2.51(2) 
(black dashed line) for random percolation. Simulations have been performed on graphs with 
= 8192 000. 




Figure 5. Explosive percolation on random SF networks. Percolation threshold Pc{N) as a 
function of the degree exponent A for different network sizes A^. The black line is the infinite 
size limit extrapolation of the critical threshold, performed by applying Eq. ([6]). The percolation 
threshold vanishes for A < Ac ~ 2.3, while it is non-zero for A > Ac. 
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Figure 6. Explosive percolation on SF networks for degree exponents A = 2.5 (top diagrams) 
and A = 2.8 (bottom diagrams). (Left) Percolation strength P as a function of the system size 
N for three different values of the occupation probability. Top: p = 0.0529 (violet diamonds), 
p = 0.0629 (orange circles) and p = 0.0729 (grey squares). Bottom: p = 0.1229 (violet 
diamonds), p = 0.1329 (orange circles) and p = 0.1349 (grey squares). The dashed lines indicate 
the best fits obtained at the critical points p = Pc = 0.0629(1) (top) and p = pc = 0.1329(1) 
(bottom), from which we get /3/i/ = 0.59(1) (top) and = 0.50(1) (bottom). (Right) The 
average cluster size S as a function of the network size A'^ for the same values of p used for P. 
The dashed lines have slopes j/u = 0.24(1) (top) and j/u = 0.42(1) (bottom). Simulations 
have been performed on graphs with N = 8 192 000. 



nodes having low degree and a small subset of nodes having high degree. The tails of the degree 
distributions are often well fitted by power laws, which lack a characteristic scale, justifying the 
name "scale-free networks" given to such systems [MlEHUMlEZlEHlEHlEOlEI]- We indicate the 
exponent of the power law degree distribution with A. Nodes with large degree, called "hubs", 
have a fundamental role for the structure and dynamics of networks. For instance, they "keep" 
a large portion of the system in the same connected component, guaranteeing stability and 
compactness to the structure. In particular, in random scale-free (SF) networks with degree 
exponent A < 3, there are many hubs and a very small fraction of links (vanishing in the limit of 
systems of infinite size) manages to keep a macroscopic fraction of nodes of the graph together 
in the same connected component. This amounts to saying that the percolation threshold is 
zero laElElElE]. 

The Achlioptas process with product rule for SF networks is implemented as follows. We 
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Figure 7. Explosive percolation on SF networks for degree exponent A = 3.5. (Left) Percolation 
strength P as a function of the system size for three different values of the occupation 
probability: p = 0.2214 (violet diamonds), p = 0.2224 (orange circles) and p = 0.2234 (grey 
squares). The dashed line indicates the best fit obtained at the critical point p = Pc = 0.2224(2), 
from which we get /S/v = —0.06(3). (Center) The average cluster size 5 as a function of the 
network size for the same values of p used for P. The dashed line has slope j/v = 0.40(9). 
(Right) Cluster size distributions Ug measured at the critical point for PR (grey squares) and 
random percolation (orange circles). The exponents are r = 2.2(1) (red dotted line) for PR and 
r = 2.94(1) (black dashed line) for random percolation. Simulations have been performed on 
graphs with = 8 192 000. 



start with a set of N nodes and a given degree sequence {ki, k2, ■ ■ ■ , k]\f}. The degrees of the 
sequence are extracted from a power law distribution with exponent A. The average degree (A;) is 
set equal to 5. Every node i holds ki stubs (half-links), where ki is the degree of i. By attaching 
such stubs in pairs one builds a network with the desired power law degree distribution with 
exponent A. If stubs are randomly connected we would have the simple configuration model [32J, 
which yields a random SF network. For the Achlioptas process at each iteration two pairs of 
stubs are selected and the PR determines which pair of stubs has to be eventually joined in a 
link (the PR applies as shown in Fig. [1]). 

The first remarkable result concerns the threshold. For random percolation it is zero for 
A < 3 and non-zero for A > 3 [5]. However, for explosive percolation, the threshold is non-zero 
already for A > A^ ~ 2.3 (Fig. [5]) [121 [l3] . 

In Fig. [6] we show the scalings at pc of P and S for A = 2.5 and 2.8. Here we see an important 
difference with the cases of lattices and Erdos-Renyi random graphs, namely that the scaling of 
P at Pc is non-trivial, as P decreases with A^ as a power law in both cases. This is very different 
from what one would find for a discontinuous transition, in which P would be approximately 
constant with A^. For S we also find power law scaling (Fig. [6l right panels), as we had previously 
seen for lattices and Erdos-Renyi random graphs. Overall, the transition looks like a standard 
continuous transition. 

For A > 3, however, the scenario changes. In Fig. [7] we show the results of our finite size 
scaling analysis for SF networks with exponent A = 3.5. For random percolation on SF networks 
it is well known [33] that (3/i^ = 1/(A - 1) and 7/1/ = (A - 3)/(A - 1) for 3 < A < 4. For A > 4, 
we are in the mean field limit and the exponents are independent of the degree exponent A: 
f3/i/ = ^jv = 1/3. These values coincide with the exponents for the percolation transition of 
Erdos-Renyi random networks. SF networks are very similar to Erdos-Renyi random networks 
in the limit A — t- 00. The explosive percolation transition looks as on lattices and ER random 
networks. The scaling of P at p^ is trivial (Fig. [3 left), with (i/v = —0.06(3), which is essentially 
zero. The power law scaling of S at pc is non-trivial (Fig. [71 center), with 7/1/ = 0.40(9). The 
Fisher exponent r = 2.2(1) (Fig. [71 right). 
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Table 1. Percolation thresholds and critical exponents for random percolation (RP) and 
Achlioptas process with product rule (PR). 



4. Summary 

We have found two different classes of phase transitions for explosive percolation, that we shall 
call class A and B. Class A includes ER random networks (the original system studied by 
Achlioptas et al. [TO]), lattices and SF networks with degree exponent A > 3. Here we have 
observed an apparent saturation of the order parameter P at pc with the size of the system N , 
due to an exponent (3 that is very small, compatible with zero. Class B includes SF networks 
with A < 3, where there seems to be no room for a discontinuous transition, as the results of 
the finite size scaling analysis are fully compatible with a standard continuous transition. Still, 
apart from the anomalous scaling behavior of the order parameter P at pc observed in class 
A, we found that the other percolation variables display power law scaling in all cases, without 
exceptions, just as in second-order phase transitions. The most striking feature is the fact that 
the size distribution of the "finite" clusters at pc is a power law, not exponential or Gaussian 
as one expects from first-order phase transitions. 

Da Costa et al. [22] have recently proved that an aggregation process in the same spirit as 
Achlioptas processes leads to a continuous phase transition, characterized by an exponent j3 for 
the order parameter P that is very small, (3 ~ 0.0555. They also demonstrated that the fact 
that the transition is continuous for the special process they considered enforces the same type 
of transition for Achlioptas processes too, with generally diff'erent but still small values for (3. 
This is fully consistent with what we have found for the systems of class A. Da Costa et al. have 
derived important relations for the critical exponents. For instance, (3/u = /3/(4/3 + 1) and the 
Fisher exponent r = 2 + /3/(3/3 + 1). For /3 ~ 0.0555 the exponents' ratio P/v ~ 0.0455. This 
value is compatible, within errors, with the values of (i/v computed for the systems of class A 
(see Table 1). So, the apparent saturation of the order parameter P with the system size N could 
be indeed due to the smallness of the exponent /3. Moreover, the Fisher exponent r for every 
transition we have investigated is very close to 2. From the relation r = 2 + /3/(3/3 + 1) < 2 + /3, 
and the fact that (3 is small, we see that r will always be very close to 2, just as we found. The 
difference between r and 2 is hard to determine numerically. Likewise, the values of r derived 
from different Achlioptas processes can be hardly distinguished from each other. We remark 
that this holds for lattices, Erdos-Renyi graphs and SF networks with degree exponent A > 3, 
while the explosive percolation transition on SF networks with A < 3 appears as something very 
different (/3 is not small here), and deserves further investigations in the future. 



We have also found that there is a non-zero explosive percolation threshold for SF networks 
for A > Ac ~ 2.3 (Fig. [5]), while the threshold for random percolation is zero as long as A < 3. Cho 
et al. [12j suggested that this is due to the non-random addition of links during the Achlioptas 
process, because of which the degree distribution of the system during the growth deviates from 
that imposed by construction, which will be eventually reached at the end of the process. This 
has been indeed verified numerically |12^ [T5] . 
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